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equations. The method is demonstrated using examples from 
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QUASI -GENERALIZED VARIABLES 
J. Baumgarte, G. P. Ostermeyer 

1. Introduction / All * 

Systems of mass points and fixed bodies, which are subject 
to certain secondary conditions, occur in analytical mechanics. 
These are idealized descriptions of real physical combinations 
between single bodies such as poles, joints, skids, cutting 
edges, etc. 

One may proceed as follows with the mathematical 
description of the system. After selecting suitable coordinates 
q , q t = (q^ , ... , q n ) , which may implicitly involve im- 

part of the combinations, the Lagrange function L is determined: 

l(q> 0 — + 'i\ ~~ f 

= t t) q+ B(q, t)q + T 0 (q, t) - T'(f/, t) (1.1) 

as are forces Q, which do not follow from L (non-potential 
forces) 

Q =Q(q>q> t) > (Oi, ... , Qn ) . d.2) 

The combinations which have not yet been considered in the 
chosen coordinates q are the so-called holonomic equations: 

/,(<*,<) =°. v = l — r (1*3) 

and the linear non-holonomic combinations of the velocities q 

9x(q, q, t) = 0 » k = \, ... ,r <n — s . (1.4) 


♦Numbers in the margin indicate pagination in the foreign text. 
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The resulting equations of motion of the observed system 
are now called 


(Lq) Lq — 0 — /gr/Jj + gqfl 2 . ( 1 • 5 ) 

Statements (1.3) and (1.4) for the coordinates, and the 
velocities, respectively, are contained herein in explicit form 
as so-called limits. 

The Lagrangian multipliers and yu x have to be calculated 
from secondary conditions (1.3) and (1.4). Normally, one 
proceeds as follows: holonomic combinations (1.3) are derived 

twice, non-holonomic combinations (1.4) are derived once totally 
in accordance with an independent variable, time t. 

h = 0 > 9y = 0 • (1 . 6 a, b) 

These relationships are then linear in accelerations q. 

If one inserts the similarly linear equations of motion in 
(1.6), one obtains an algebraic system of equations with 
Lagrange multipliers. Because of the usual independence of the 
several combinations, these multipliers can be clearly 
calculated as functions of q, q, and t. 

Although equations of motion (1.5) are very convenient in 
considering additional combinations (1.3) and (1.4), they have a 
major disadvantage with regard to numerical integration. The 
numerical solution of (1.5) exhibits characteristic errors of 
the combinations (1.3) and (1.4), which, practically independent 
of the quality of the program of integration, quickly render the 
solution useless. Essentially the errors in the combinations 
increase linearly--it should be noted that this implies a 
quadratically increasing error for coordinates q--and 
quadratically in the combinations (1.3) . This pattern of errors 
is implicit in the calculation of the Lagrangian multipliers 
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using (1.6 a, b) . Additional instabilities are transofered to 
the equations of motion. 

Although one may achieve very good results by using 
stabilization strategies (see section 2) to diminish these 
instabilities, it is actually unsatisfactory from the analytical 
point of view, out of explicit consideration for combinations 
such as those in movement equations (1.5), to transfer 
instabilities into the motion equations only to eliminate them 
later in the numerical solution. 

Theoretically these difficulties can be completely avoided, 
when all holonomic combinations are implicitly considered via 
suitably chosen generalized coordinates *q, and when the /472 

non--holonomic combinations are eliminated through "non-holonomic 
velocity parameters" (Hamel [1]). While non-holonomic 
combinations are eliminated rather easily- -non-holonomic 
secondary conditions are generally linear in velocity--this 
procedure is not always very easy or even possible for holonomic 
combinations, which are generally non-linear in coordinates. 
Therefore, the following is recommended in related literature. 
Holonomic combinations are completely derived once from time. 

Then they can be treated formally like non-holonomic 
combinations. The reduction using non-holonomic velocity 
parameters to, to r = (toj, , <*>»_,_«) then yields a system of 
differential equations of the order 2n - r - s 

q — F(q, to, t) , to — G(q, ta,t) . ( 1 . 9 a , b ) 

In the process, the derived holonomic combinations (1.3)' 
are thus implicitly stated using the variables id . The 
combinations (1.3) in the non-derived form are now first 
integrals for the n equations (1.9 a). 
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Peculiarly, the numerical solution disturbs first integrals 
as much as the explicitly considered non-holonomic 
combinations. Here too, it is true that instabilities occur 
which make themselves numerically noticeable in strongly 
increasing errors with reference to the combinations. These 
instabilities are caused by the incomplete-- though 
mathematically possible--reduction of the order of differential 
motion equations. 

It is the goal of this work to formulate motion equations 
so that, even with explicitly considered combinations, not only 
the analytical, but also the numerical, solution is error-free 
regarding the secondary conditions. 


2. Stabilization 

The point of departure for stabilizations of combinations 
with motion equations of the order (1.5) or (1.9) involves 
(Baumbarte [2]) calculating the Lagrangian multipliers in (1.5) 
no longer using (1.6 a) or (1.6 b) , but rather the 
asymptotically stable differential equations 

/; + 2«/; + P 2 f, = 0 , fi = « > 0; g x + yff x = 0 , y> 0. (2.1 a; b) 

The Lagrangian multipliers and^u^ calculated thusly, now 
contain the classic portion together with expressions which 
disappear identically in the analytical solution, but for which 
regulator terms constitute the numerical solution regarding the 
combinations. The stabilizations (2.1 a, b) limit the errors 
regarding the combinations (1.3) and (1.4) in the numerical 
solution. From the regulat ing-technical point of view, (2.1 b) 
is a P- and (2.1 a) a PD-regulator . 

Extremely precise stabilizations (Ostermeyer [3]) may be 
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obtained by inducing an I-member into the regulator, since then, 
on one hand, the temporal median of the numerical combination 
error disappears and, on the other hand, the coefficients in the 
stabilizations can be directly indicated after inducing such a 
substitute constant T* using the optimum criterion employed in 
regulator technology. 

t t 

/drj=0; g + ^^ g +J^Jgd rj=0, T*>0. (2.2a; b) 

These stabilizations certainly increase the order of a 
differential equations system--with the additional differential 
equation 

a = f(x, t) or = g(x, x,t) 

one obtains the integrations necessary with z using the 
combinations, but the substitute time-constant "t*" only 
accounts for the process of the "real", i.e., error-ridden, 
integration. 

An optimum value can be given even for quantity T* , 
whereby essentially the salient frequencies of the then adaptive 
regulator, as a function of the highest occurring frequency in 
the mechanical system, are conducted (OSTERMEYER [4]). Thus the 
parameters in the stabilizations are completely determined. 

Even first integrals can be stabilized with the integration 
of a given mechanical system (BAUMGARTE [5]). The integrals 
I,(x, x, t) = 0 are hereby interpreted as non-linear non-holonomic 
combinations, just as in the velocities. From the Gaussian 
principle of the smallest limit, one can formally obtain the 
motion equation with the "limit" 


2 r ]'* 


4 / + J 


/ + 

’* ' ^ 1C 2 


T .,j 



The usual calculation of the Lagrangian multiplicator 


from 


i = o < 2 - 3 4 > 

yields, of course, S Ot since, according to the definition, 
the motion equations 1 ~ are identically realized. However, 

r 

if the multiplier d is calculated from J v in the stabilized form 
(2.1 b) , (2.2 b) , one obtains an algebraic expression for d. 

This expression is a pure regulator-term and disappears 
identically upon the exact solution. 

The stabilizations are a simple aid for controlling 
numerical instabilities with great accuracy and little effort. 

But some questions still remain (see [4]). Although the 
procedure of stabilization seems very plausible, it has not yet 
been determined in a general fashion whether and in what sense 
the numerical solution of (1.5) with stabilized combinations is 
altogether better than the numerical solution of (1.5) without 
stabilization, if one disregards the component of the total 
numerical error, which makes itself noticeable in the 
disturbance of the combinations. 


3. Quasi-Generalized Variables /473 

A differential equation system is given in the variables 
X ~ X (t) with X T «=(*,. ... ,*„): 

<b=I ?(*,<), «€/ = [< 0 ,y, with x(t —t 0 )=x 0 eR n . (3.1), (3.2) 

Apply an additional combination to this system 

/(*,<)= 0, with /(* 0 ,* 0 ) = 0. (3.3) 
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There are two classic means in mechanics to obtain a 
solution to this problem. 

A. Reduction 

Here the generalized variable 2/(0 > U T ~ (ih Vn-i) , generalized 

on n - 1, which covers the subspace given by (3. 3) , is 
transformed. With 

x*=x(y,t), x v -f a = o (3.4a) 

one obtains the differential equations 

Xyii + t),t) , or ij = {Xy ■ .rj)- 1 x y (F — r-i) . (3.4b) 

B. Lagrangian Multiplier Technique 

Here, the minimum of the quadratic form Z 

Z =i(i -*Y(x -F) (3.5) 

with regard to the quantity 'X is sought under the secondary 
condition (3.3). Once (3.3) has been totally derived according 
to time, then one obtains a linear combination for dc. With the 
Lagrangian multiplier yU, one obtains from 

~(Z -/</)= 0 (3.6) 

the differential equation system 

i=F+fif m , /(as, <)=0. (3.7 a, b) 

# 

If one inserts (3.7a) into f = 0, one has an algebraic 
equation for calculating jjl. The integration of the n 
differential equations 
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(3.8) 


x 


= i<' 


/xF + ft 

il-U U 


yields the solution to the problem. 


Procedure A yields numerical solutions which conform 
exactly to the combinations. The disadvantage of this procedure 
is that, under certain circumstances, the search for generalized 
variables is not very easy. 

Procedure B has the advantage that one simply obtains the 
resulting differential equation system. The numeric instability 
of the combination is disadvantageous. 

Principle procedure C, submitted here--introduction of 
quasi-generalized variables, see below — no longer considers a 
differential equation system with the actual involved variables 
x and y, but rather introduces auxiliary variables x, x T = (x lt , 

which are unrelated, but which contain the involved related 
variables x (i.e. x with f(x,t) ^=0) in a simple fashion. The 
differential equation systems pertaining to ft should be easily 
constructed and formed in such a way that, after numerically 
solving the differential equations in x, one may transform from 
(3.8) into the numerical error-free solution of x regarding the 
combination. 


The point of departure is the formulation of the problem in 
the form (3.7) or (3.8). If X =X{t) is an arbitrary 
(continuous) function of time, it is valid that: 

* = F + Xf x + (ft — X) f x • (3.9) 

With the transformation 

x^x-XU ( 3 . 10 ) 
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one obtains further 


tc=F-X(f x y +{(i-k)f x . (3.11) 

If one executes the differentiation of the terra (f ) and 

X 

re-inserts (3.7), (3.11) is simplified to 

x ~ (i — ;./**) (f + fiU) - A/* (3.12) 

Since on this premise, the analytical solution x of (3.7) 
or (3.9) of the combination conforms exactly to f(x,t) = 0, the 
following is valid using (3.10) 

ftx, t) = f(x(x, A), t) = 0 . (3.13) 

Using (3.10), the right side of (3.12) can be completely 
formulated in the new variable X. According to the 
construction, x for 1 ^ 0 no longer lies on the hypersurface 
given by the combination. 

Until now A can be set arbitrarily. The solution of /474 

(3.12) in the new variable X contains the first integral 

(3.13) £ [X)\ t,) ~ O' How should A be suitably chosen? It is 
possible, here for instance, to set A 0. Then, * = .x-fc • /* 
with the set constant c is valid for the analytical solution. 

If one is to insert the constant c for X into the equation 
resulting from (3.10) 

X = (7 — ;./«) (F + Fix) -lf*t (3.14) 

the calculated solution for x obtained from the numerical 
solution for X in (3.14) using (3.10) disturbs the combination 
as much as the numerical solution of the starting system (3.7) 
or (3.8). Up to now one has only executed a coordinate 
transformation which leaves the same numerical difficulties as 
before, because the numerical solution of (3.14) itself disturbs 
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its first integral (3.13). But only the exact satisfaction of 
(3.13) with the integration of (3.14) guarantees that the 
re-transformation of arbitrary values for 5 : to the. starting 
coordinates x excludes errors in the combination f(x,t) = 0 or x. 

The first integral (3.13) can be used to eliminate and 
in (3.12). This reduction (of the arbitrary choice of X) has 
the result that (3.13) is implicitly considered in (3.12) and 
that: thus (3.13) can no longer be disturbed by the numerical 
solution of (3.12). Furthermore, using simple calculation, one 
finds the identity 

X zz- jx y (3.15) 

so that the motion equations now have the following form 

*==(/- 1/**) F- Wf'mm •/*+/*(). a =A(5,0 (3.16 a, b) 

# 

The analytical expression for X can now be calculated from the 
total derivation of (3.16 b) to t, or, for instance, from the 
total derivation of (3.13) to t. (This is equivalent to the 
transformation of jx in (3.7), (3.8) into x,X ) • The first method 
presumes that X can be indicated explicitly from the implicit 
system of equations (3.10), (3.13). The second method yields an 
explicit expression for X even when X is numerically 
iteratively established (see sections 4, 5) . 

This choice of procedure for X corresponds in principle to 
the possibility to partially integrate the Lagrangian 
multiplier jj in (3.7a). (See Baumgarte, Ostermeyer [6]). 

In this sense, n variables X are quasi-generalized 
variables, since the differential equations (3.16) formulated by 
them no longer contain the classic "limit" term juj x and since 
they implicitly contain the relevant solution x, which exactly 
satisfies the combination after the construction. 
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The last point is given inasmuch as the first integral 
(3.13) regarding the variable^ through the reduction regarding 
the freedom in the choice of X is also exactly numerically 
satisfying . 

C. Procedure for Induction of Quasi-Generalized Variables 

The new variable x, x T = (x v ... , *„) on n is here transformed 

with 

x — x - A/ x , f(x, t) = 0 , x - x(x, X, t) , X = X(x, t) . 

The differential equation system is 
x = (I - X) xx ) F - XQUU + fxt) ■ 

These are n first-order differential equations in X. 
relevant solution x is established after the (numerical) 
integration of (3.18) , x is error-free regarding the 
combination f(x,t). The starting quantities of x are calculated 
from Xq using (3.17), and calculated using an arbitrary 
value X(t 0 ) . For the choice Xfto) —0 it is true that 

x(t = t 0 ) = x(t = fo) . (3.19) 

4. Treatment of First Integrals 

If, for a differential equation system 

x=--F(x,t) (4.1) 

the combination f(x,t) = 0 is a first integral (f is then a 
so-called inner connection), then fjc +} t ~ f z F + /, = 0 is valid, 
and the Lagrangian multiplier is 0 with the tranf ormation 
into quasi-generalized variables. The differential equation 
system (4.1) is transformed with x=x+Xf x into 


(3.17) 

(3.18) 
The 
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X =■ (I - Xj xx ) F - Xf xt 


(4.2) 


Example 1: The Quaternion Differential Equation 

While describing the position of a fixed body or a 
gyroscope, with a given turning vector 

<0=0 >(t) . Oi T — (ftjj, to 2 , <0 3 ) 

singularities occur in Euler's angles ("frame-locking"). In 
order to avoid these difficulties, other position coordinates 
are often introduced. 

The turning quaternions a, x T = fo, x 2 , x 3 , z 4 ). constitute /475 

one possibility. The related differential equation is . t , 

x — \ Six ( 4 . 3 ) 

with 

( 0 —w 1 -w 2 

Wl 0 
W 2 ~ m 3 0 

0> 3 (0;j — CU X 

The quaternions satisfy the first integral 

/(*)=**- 1 =0. (4.5) 

Numerical integration of the system (4.3) disturbs the 
combination (4.5) linearly in t. The procedure C (see section 
3) leads, using the transformation 

x == x - ).f x = x — 2)jx , x = 5/(1 - 21) (4.6) 

to the new differential equation system 
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(4.7) 


x = i {I -2X1) Q- (1 ~ 32 X[ bzw ‘ ' * = v Sis5 ■ 

The differential equation system in the quasi-generalized 
variables (4.7) has the same form as the starting system. 

For the re-transformation into the starting variable, X is 
needed. The insertion of (4.6) into (4.5) yields: 

A = ->-(1-1*1). (4.8) 

Therefore 

x = a$/|e| • (4.9) 

The rule used to divide quaternions by their norm, after 
every step of integration, to take the next step of integration 
with the new starting values for the quaternions x, has been 
known for a long time as a numerical variant of the integration 
of the differential equation (4.3). (4.9), however, signifies 

more: when dealing with the problem (4.3), formulated in the 
quasi-generalized variables (4.7), one may not interfere in the 
integration. One only has to realize the depiction onto the 
starting quaternions at those points of time when the solution 
is relevant. The resulting x vector is, independent of these 
points of time, a numerically correct solution of (4.3), which 
exactly satisfies its inner combinations. 

Example 2: The Unperturbed Kepler Problem 

With x = {x l> x v x i ) and r := jxj , the motion equations are 
K s 

i-f — jc=0, k* = parameter of gravitation (4.10) 

The analytical solution to the unperturbed Kepler problem 
can be given immediately. In practice, one will often work with 
a mildly perturbed Kepler problem (see Stiefel, Scheifele [8]). 
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However, the numerical difficulties with integration of Kepler 
equations may be examined with particular ease in the 
unperturbed problem. 


An essential difficulty with the integration of (4.10) is 
the upset of the energy integral 


/ 




E = const. < 0 . 


(4.11) 


Since the energy determines the frequency of revolution (j, 
Kepler's orbit will be only mildly perturbed with numerical 
errors in (4.11)--the Kepler problem is stable in orbit--but the 
error in calculating the location of the satellite in orbit 
increases quadratically with time. 

The consideration of (4.11) using quasi-generalizeed 
variables in procedure C can be executed when (4.10) is written 
as a first-order differential equation system. 



(4.12 a, b) 


with 


/(*,»)=■ 



(4.12 c) 


The transformation could be executed with regard to x and v. 
Consequently, one has a form of the gradient procedure. 


However, the principles of mechanics indicate that-- (4.11) 
being understood as a non-holonomic combinat ion--the combination 
is only considered in the equations (4.12 b) , meaning that only 
the quantities v have to be transformed into the quasi- 
generalized variable xr. 


Using 


/4 76 


/ v '= v — Ac , v = »/( 1 — A) 


(4.13) 


16 



can be calculated from (4.12 c) , using the abbreviation V for 
the potential 

F = - K*l\x\ (4.14) 

the multiplier A into 


A = 1 - \v\t\'2(E - V) 

and thus the transformed motion equations into 


(4.15) 


X - 




M 


A 2 


\'2{E - V ) r 


■ x . 


(4.16 a, b) 


Re-transformation into the variables (x, v) yields 
solutions to (4.12) which exactly satisfy the energy integral. 


Example 3: The Perturbed Kepler Problem 


Generally, one always has to work with a mildly perturbed 
problem: 

A 4 

X x = EQ{x,x,t) . (4.17) 

If the forces are potential, the energy integral (4.12 
c) , with additional potential terms, continues to be valid. The 
transformed equations keep the structure (4.16). If the 
terms £ Q are non-potential forces, the energy integral is no 
longer valid. With expansion, however, of the motion equations, 
the transformation can even be executed here. Using 

x — n , b = —V m + iU, E =eQv (4.18 a, b, c) 


one again induces a first integral: 



As the errors in the additional scalar differential 
equation (4.18 c) are significantly smaller, using the numercial 
interpretation, than the errors in the system (4.18) with regard 
to the integral (4.18d), the transformation into quasi- 
generalized variables yields, even here after re-transformation, 
significantly better results. The transformed system reads 




M 


)/2(E - V) 


( - V x + cO) , E = e(Q f2(E - F) . 


(4.19) 


Example 4: The Symmetric Gyroscope 


With the principle moments of inertia = A,, 1 3 , which 

belong to the center of gravity (fulcrum) , and with the 
components of the angular velocity with regard to the solid 
principle axis of inertia, Euler's equations read 

JTjG) t = Koj^ j A 1^2 “ — JSo ) i j / 3 CO 3 — 9 to 3 — fOg , A = j cOg . (4.20 a, b, c) 

■*3 

The system contains the first integral 
/(to) == Y /,(«? + w|) — E — 0 . 


Using the transformation 

COj °>2 

Wi ' _ (l -U~) 

follows, from (4.20), the system 

IjWj == h<0 2 , I ^( 1)2 “ —Kg)} 1 -^3^3 = 0 . (4.21) 

Here again, one obtains the starting system, formulated 
with the new variables. The re-transf ormation reads 



i 


2 E 

h{a>{ + to|) 


U)i 



(4.22) 
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Concerning the practically important problem of the 
perturbed gyroscope, one can proceed similarly to deliberations 
of the perturbed Kepler problem. With the vector of 
perturbation n = (N v N t , N 3 f , the system now reads 

IjWj = I\w 2 + eNy , JjWj, = — /fwj + eN 2 , K = — <u 3 , m 3 ^- const . , (4.23) 

^3 

I 3 co 3 =: -{-eN 3 , E — e(N 1 w l -f- N 2 co 2 ) . 

The integral is 

/ =-3-A(®l + r»5) -E = 0. (4.24) 

The transformation into quasi-generalized /All 

variables } leaves the structure of system (4.23) 
unchanged, as above. Re-transformation (4.22) into the 
values cbj, after numerical integration of (4.23), is 
error-free after construction with regard to (4.24). 


5. Concerning Numerics with Non-Explicitly Realizable 
Transformation into Quasi-Generalized Variables 

In all the previously examined examples, the transformation 
could be explicitly executed. The non-symmetric gyroscope 
constitutes an example with which this no longer seems 
possible. For purposes of simplification, let us examine it 
here in its unperturbed form. The perturbed case can easily be 
examined as demonstrated above. 

Example 5: The Non-Symmetrical Gyroscope 

With the terms of the previous example, Euler's equations 
here read 


— (1* — I») > / s a> 2 = (J 3 — - / x ) co 3 oj t , I 3 6) 3 = (/j — I 2 ) (OiCo., . (4.22 a, b, c) 
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The system contains two first integrals: 

A = t<A«>! + Vf + M -<A = 0, h=\ {l\w\ + l\ m \ + !§«>f) - C 8 = 0 . (4.23), (4.24) 

Here and C 2 are constants for the unperturbed case. 

The specification for transformation reads 


( 0* — (Of — XJiWi — A 2 /fa)i , i -■ 1, 2, 3 . (4.25) 

The motion equation with quasi-generalized variables reads 


(I — — j-gjjl) (A A) ^jj^3 

(I-V2-VD U-V3-- Vi)‘ 


(4.26) 


The remaining equations can be obtained via cyclical 
permutations. 

In (4.26), however, the Lagrangian multipliers , A a are 
still to be calculated using (4.23) and (4.24). This 
calculation can be performed numerically. One must consider, 
however, that with every new call-up of the differential 
equation (4.26), through the numerical routine of integration, 
the actual multipliers 

A j.i A) — 0 , Utfi* it A> A) — 0 > A — A(®<) » A = A( w *) (4.27) 


should be established via an iteration procedure. 


It may also occur that not only the Lagrangian multipliers 
X, but even the initial variables x as a function of the 
quasi-generalized variables X cannot be explicitly established. 
In this case, the (n + m) equations 


x — x — f n \ — 0 , f{x, t) —0 , / a =- (A, ••• 1 1m) • — (A> ••• > A») (4.28 a, b) 

have to be solved iteratively into x and X to establish the 
differential equation system 
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for the integration routine,. 


To avoid difficulties while searching for the zero digit, 
as well as long iteration periods, the following procedure has 
proved effective. 

The starting quantities of the quasi-generalized variables 
are defined with only the exception of a free choice of starting 
quantities of /l. For purposes of simplification, A (t—t 0 )—o 
is given (see (3.19)). If the transformations into x and X 
following from (4.28) can be explicitly executed, it will be 
evident that the Lagrangian multipliers increase, in general, 
linearly when using the numerical integration of (4.29), and 
even quadratically , when using the quasi-generalized variables 
for holonomic combinations (see section 6) . 

If one now evaluates the system (4.28) during and after the 
first step of integration At, the desired solution vector x(t 0 + At) 
is present. At the beginning of the next step of integration, 
the vector value X can again be designated zero, A(< 0 + At) = 0. 
With the vector of solution »(<„ -f At 0 ) , one obtains, using 
(4.28), the initial quantities x(t 0 + At) for the next step of 
integration. 

If one proceeds in such a way through every step of 

integration, vector X remains very small. Numerical experiments 

have indicated that, using this procedure to evaluate equations 

(4.28), only one or two steps of a Newtonian procedure are 

needed to obtain very high accuracies--for instance, to obtain 
5 10 

the factor 10 or 10 times better than the stabilization 
in section 2. Using our integrational procedure, RKF 7 [8], the 
differential equation system is evaluated at a total of 13 times 
per step of integration. Having executed each time the 
iteration of (4.28), the resulting calculation time was slightly 



longer, by a factor of 1.5 than the integration of the starting 
differential equation system. 

6. Treatment of Holonomic Connections 

The holonomic combination is f(x,t) = 0. If one examines 
it in derived form like a non-holonomic combination using 
non--holonomic velocity parameters, one obtains the system 
(1.9). f(x,t) = 0 is then a first integral of equations (1.9) 
and can be treated with the procedures in sections 4 and 5. 

If f(x,t) = 0 is examined using a limit in the /478 

motion equations (see (1.5)), these read;: 

M(x, 1) x = F(x, x, t) + nf Xr 
or, as a system of the first order 

i = F(x, vt) + n) x . (6.1 a, b) 

In this instance, transformation into a quasi-generalized 
variable is also possible. Using 

x: = x —XU (6.2) 

one again obtains, using the secondary condition 

x = x(x, X, t) , X = X(x, t) . (6.3a, b) 

In Baumgarte, Ostermeyer [7] it is shown that conversion 
into quasi-generalized variables X converts the equation i6.1) 
into the following system 

* = M~'(v + X{fxY) , ® = F + W*)' (6.4 a , b) 

usina c. = ® + Xf x , (/*)' = “ Vx) + /**• 
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It is given that AiEEyU. In this instance, the function x, 
established via numerical solution of (6.4) and (6.3), is also a 
solution to initial system (6.1) which is error-free with regard 
to the combinations f(x,t) == 0. In the system (6.4), A can be 
calculated either by total derivation of A or from the derived 
combination f = 0. Since the derived combination is linear in 
its velocities, transformations into v can always be explicitly 
executed; 

«.= «(«, K <M) = v - \f x , i = A(», x, t) ~ tfEjLh . . 

Jz ’ /* 

If this is substituted into (6.4), only the calculation of 
X 6c, Xt) and A = \(x,t) remains. If this is not explicitly 
possible, one may proceed numerically in this instance, as 
demonstrated in section 5. 
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